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ABSTRACT 

Context. Thermal wind emission in the form of free-free and free-bound emission is known to show up in the infrared and radio contin- 
uum of hot and massive stars. For OB supergiants with moderate mass loss rates and a wind velocity distribution with /? ^ 0.8 . . . 1 .0, 
no influence of the wind to the optical continuum, i.e. for A S l.Oyum, is expected. Investigations of stellar and wind parameters of 
OB supergiants over the last few years suggest, however, that for many objects fj is much higher than 1.0, reaching values up to 3.5. 
Aims. We investigate the influence of the free-free and free-bound emission on the emerging radiation, especially at optical wave- 
lengths, from OB supergiants having wind velocity distributions with p > 1.0. 

Methods. For the case of a spherically symmetric, isothermal wind in local thermodynamical equilibrium (LTE) we calculate the free- 
free and free-bound processes and the emerging wind and total continuum spectra. We localize the generation region of the optical 
wind continuum and especially focus on the influence of a yS-type wind velocity distribution with /3 > 1 on the formation of the wind 
continuum at optical wavelengths. 

Results. The optical wind continuum is found to be generated within about 2 R, which is exactly the wind region where /? strongly 
influences the density distribution. We find that for /? > 1, the continuum of a typical OB supergiant can indeed be contaminated 
with thermal wind emission, even at optical wavelengths. The strong increase in the optical wind emission is dominantly produced 
by free-bound processes. 

Key words. Stars: early-type - supergiants - Stars: winds, outflows ~ Stars: mass-loss - circumstellar matter 



, 1. Introduction 

It is well established that for massive stars the appearance of 
(thermal) excess emission at infrared (IR) and radio wavelengths 
■ is caused by free-free and (to a small fraction also by) free-bound 
I emission generated in their winds (see e.g. Panagia & Felli l 19751 
Olnon'1975). At radio wavelengths, the free-free excess emis- 
sion is usually used to derive the mass loss rates of hot and mas- 
sive stars (see e.g. Lamers & Leitherer |1993t Puis et al. |1996l ). 

Waters & Lamers (1984) have investigated this excess emis- 
sion in detail, especially in the near-IR region. These authors 
studied the influence of free-free and free-bound emission to the 
total continuum and pointed already to the importance of the 
wind velocity (and hence the wind density) distribution that can 
severely alter the wind contribution in the IR. 

However, the investigations of Waters & Lamers ( 1 19841 1 were 
restricted to the IR and radio range, i.e. they calculated the wind 
contribution for A k 1 fj.m only, while during earlier studies 
Brussaard & van de Hulst (1962) had noted that free-bound pro- 
cesses might become very important in the optical and UV range 
for temperatures typically found in the winds of hot stars and su- 
pergiants. 

Whether the free-bound emission in the wind indeed influ- 
ences the optical continuum, depends severely on the density 
distribution. For line-driven winds, the density follows from the 
equation of mass continuity, i.e. it is proportional to the mass 
loss rate, and inversely proportional to the wind velocity. A high 



wind density can therefore be reached by either a high mass loss 
rate, or a low wind velocity. 

For massive stars with pronounced high mass loss rates like 
Wolf-Rayet stars. Luminous Blue Variables, or the group of B[e] 
stars, it is well known that the wind not only influences, but even 
dominates, the optical spectrum. In these objects, the wind is 
usually optically thick even in the visual range and thus com- 
pletely hides the stellar spectrum. Some recent examples of the 
thermal wind influence in the form of free-free and free-bound 
emission at optical wavelengths have been published, e.g. by 
Guo & Li (2007) for the case of Luminous Blue Variables, and 
by Kraus et al. ( I2007I I for a Magellanic Cloud B[e] supergiant. 

The second density triggering parameter is the wind velocity 
distribution. For OB-type stars with line-driven winds, the wind 
velocity is very often approximated by a so-called yS-law, where 
/3 is in the range of 0.8 ... 1 .0. Such a /3 value causes a rather fast 
wind acceleration at the base of the wind, and the wind reaches 
its terminal value within a few stellar radii (see e.g. Lamers & 
Cassinelli |1999l ). Consequently, the region of very high density 
that might cause enhanced free-bound emission is restricted to 
an extremely small volume around the stellar surface. We can, 
therefore, expect that for OB-type stars, which have only moder- 
ate mass loss rates and whose wind velocity distributions have /3 
values between 0.8 and 1.0, there will be no noticable influence 
of the wind on their optical continuum emission. 

During the last few years, huge effort has been made to de- 
termine precisely the stellar and wind parameters of OB-type 
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supergiants in the Galaxy (e.g. Kudritzki et al. 119991 Markova 
et al.'SOOT; Fullerton et al. '2006'; Crowther et al. |2006l Puis et 
al. 12006 ), in the Magellanic Clouds (s ee e.g Evans et al. 120041 
Trundle et al. 2004; Trundle & Lennon |20051 l. and beyond, e.g., 
in M31 (e.g. Bresolin et al. 120021. Interestingly, many OB su- 
pergiants are found to have rather high /? values (up to fi - 3.5, 
see Sect. 13. 21 ). Such high y6-values, resulting in a strong density 
increase close to the stellar surface due to a much slower wind 
acceleration, should have noticable effects on the thermal wind 
emission not even in the near-IR, but extending also to optical 
wavelengths. In this paper, we, therefore, aim to investigate and 
discuss in detail the influence of high fi values on the wind emis- 
sion of OB supergiants, especially at optical wavelengths. 



2. Description of the model OB supergiant 

The calculation of the continuum emission of a typical OB su- 
pergiant is performed in three steps: (i) first we calculate the 
stellar emission of the supergiant with no stellar wind, (ii) then, 
we calculate the emission of the wind with the stellar param- 
eters as boundary conditions, (iii) and finally, we combine the 
two continuum sources whereby the stellar emission still has to 
pass through the absorbing wind. A justification of the use of 
this so-called core-halo approximation together with a discus- 
sion of several other assumptions and simpUfications are given 
in Sect.i] 



2.1. The stellar model 

To simulate a typical OB supergiant we adopt the following 
set of stellar parameters: Tetj = 33 000 K; - 17.2 
logL./Lo = 5.5; and \ogg - 3.4. With these parameters, we 
compute the stellar continuum emission of a hydrogen plus he- 
lium atmosphere, given by the Eddington flux, Hy. These calcu- 
lations are performed with the code of Kubat (2003 and refer- 
ences therein), which is suitable for the calculation of non-LTE 
spherically-symmetric model atmospheres in hydrostatic and ra- 
diative equilibrium. 



2.2. The wind model 

For simplification, we assume that hydrogen in the wind is fully 
ionized while further contributions to the electron density dis- 
tribution from, e.g., helium and the metals, are neglected. This 
means that for a given mass loss rate the real number density of 
free electrons is underestimated. This leads to an underestima- 
tion of the total wind emission generated via free-free and free- 
bound processes. However, since we are interested only in the 
effect on the emerging wind emission caused by different veloc- 
ity distributions, such a simplification is reasonable. We further 
neglect electron scattering, and describe the wind zone follow- 
ing Panagia & Felli (1975 ) with a spherically symmetric sta- 
tionary model in LTE. We show the frame of reference used for 
our computations in Fig.[T] The electron number density distri- 
bution, ne(r), then equals the hydrogen number density distribu- 
tion, nii(r), which follows from the equation of mass continuity. 
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This equation relates the density at any location r in the wind 
to the mass loss rate, M, of the star and the wind velocity, u{r). 
The parameters yU and niu are the mean atomic weight, for which 




Fig. 1. Frame of reference used for the calculations (after 
Panagia & Felli |1975l l. The x-axis is perpendicular to the 
drawing-plane. 



we use a value of 1.4 (i.e. solar composition), and the atomic 
hydrogen mass, respectively. 

The velocity increase in line-driven winds of hot stars is 
approximated with a y6-law of the form (see e.g. Lamers & 
Cassineni [T999] l 
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where /?, is the stellar radius, is wind terminal velocity, and 
P describes the 'steepness' of the velocity increase at the base of 
the wind. The term i>o gives the velocity at the base of the wind, 
i.e., at r = 7?,, for which we use the thermal hydrogen velocity 
given by the wind temperature. 

The wind temperature for hot stars is known to drop quickly 
with distance (see e.g.. Drew 1 1989b . and a useful description of 
the radial electron temperature distribution as a function of ef- 
fective temperature and wind velocity has been derived by Bunn 
& Drew (1992i ) 



re(r) = 0.79 Teff- 0.51 ^Teff. 



(3) 



This relation states that in the vicinity of the stellar surface, the 
electron temperature is roughly O.STeff, a value that has been 
found and confirmed by other investigators as well (see e.g. 
Drew [T989l de Koter [T993l Krticka|2P06). 

For the purpose of our investigation, it is sufficient to restrict 
the model computations to an isothermal wind, and throughout 
this paper we will use the value of - Te(R,) = 25 760 K 
(which is lower than O.STeff) following from Eq. (13) for our star 
with Teff = 33 000 K. A justification for these assumptions and a 
discussion of their reliability are given in Sect. |4.1| 

The assumption of a constant wind temperature allows for a 
simplified treatment of the radiation transfer and therefore of the 
intensity calculation. Consequently, at each impact parameter, ^ 
(see Fig.lTJ, the intensity of the wind emission is given by 



/v = Bv(re)(l-e-"'<«) 



(4) 



The optical depth is defined as the line-of-sight integral over the 
absorption coefficient of the free-free and free-bound processes. 
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with the integration limits 



and 



for 0<( <R,, 
for R,<(< Rout : 



(6) 
(7) 

(8) 



where /?out is the outer edge of the ionized wind. The total ob- 
servable flux at earth of the wind continuum emission follows 
from the integration of the wind specific radiation intensity (HJl 
over the wind zone projected to the sky 



^,wind ~ 



27T 



(9) 



where ^n,ax = ^out, and d is the distance to the object. 

With the wind absorption along each line of sight, i.e., each 
impact parameter the stellar flux passing through the wind 
zone can be calculated from. 



^ v,star ~ ^2 ^ ; 



(10) 



with the stellar intensity ly = 4-Hy. For simplicity, we restrict 
the calculation of the wind attenuation to the direction ^ = 0, 
only. This delivers a lower limit to the real attenuation so that 
we slightly overestimate the stellar flux leaving the wind zone. 

Finally, the total continuum emission is the sum of attenuated 
stellar and wind contribution, i.e. 



: ~^ -^v.wind ■ 



(11) 



2.2.1 . Gaunt factors for free-free and free-bound processes 

The absorption coefficient of the free-free and free-bound pro- 
cesses, Kv((, s) (in cm''), is given by (see e.g. Brussaard & van 
de Hulst [Tgg2] i 



Ky((, s) ^ 3.692-10' 



{l-e-^)(gg(y,T)+f(v,T))(12) 



where gff(v, T) is the Gaunt factor for free-free emission and the 
function /(v, T) contains the Gaunt factors for free-bound pro- 
cesses. To calculate the absorption coefficient, we need to spec- 
ify the Gaunt factors for both the free-free and the free-bound 
processes. 

In the literature there exist several approximations for the 
calculation of the free-free Gaunt factors in either the short or the 
long wavelength regime. To allow for an appropriate transition 
from one regime to the other (see e.g. Waters & Lamers 1984; 
Krausj2000) we use in the long-wavelength region the relation of 
Allen 0973), and in the short-wavelength region the expression 
of Gronenschild & Mewe ( I1978I I. The resulting Gaunt factors 
are in good agreement with those calculated from the approxi- 
mation given by Mihalas ( 1967), which is based on calculations 
ofBerger( |T956l ). 

In the top panel of Fig.|2l we plotted the free-free Gaunt fac- 
tors calculated over a large range of wavelengths and for dif- 
ferent electron temperatures. The curves converge to unity for 
wavelengths shorter than ~ 1 fim, while they increase with wave- 
length in the IR and radio range, where free-free processes are 
known to dominate the continuum emission. The free-free Gaunt 
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Fig. 2. Top: Gaunt factors for free-free processes. They con- 
verge to unity in the optical, but increase with wavelength in 
the IR and radio range. Middle: The function /(v, T) contain- 
ing the Gaunt factors for the free-bound processes. The values 
of /(v, T) increase steeply in the optical, while for A > I yum 
they drop quickly and become unimportant compared to the free- 
free Gaunt factors. Bottom: Absorption coefficient of the free- 
free and free-bound processes (see Eq. (fTSll). Clearly visible is 
the growing influence of the free-bound processes in the opti- 
cal range, especially for decreasing temperature. The line styles 
indicate the electron temperature. 



factors thereby depend only weakly on the chosen electron tem- 
perature. 

The situation is completely different for the function /(v, T), 
which itself is not a Gaunt factor but that contains the Gaunt 
factors of the free-bound processes. We calculate this function 
with Eq. (25) of Brussaard & van de Hulst ( 119621 ). The factors 
g„ entering this equation are thereby the Gaunt factors for the 
transitions from the free level E = hvo/n^ + hv to the bound level 
E - IivqIvP-. These individual Gaunt factors can, in principle, 
be calculated with the exact formula given by Menzel & Pekeris 
(il935l ) or by the approximation provided by Mihalas (1967). But 
it turns out, that these Gaunt factors have values between 0.8 and 
1.1 in the wavelength range of our interest (i.e. for /I ^ 0.1 ^m) 
and converge to unity for A ^ lOyum (see Fig. 8 of Brussaard 
& van de Hulst |T962l l. It is therefore reasonable to use g„ ^ 1 
for all levels of interest. With this approximation, we calculated 
the function /(v, T) over the same wavelength range and for the 
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Tj„ = ,13 000 K 
R. = 17.2 Rg 
M = 5-10"'^ Mg/yr 
V, = 1400 km/s 
(3 = 1 .0 



Wavelength A 

Fig. 3. Continuum emission of a typical OB supergiant (solid 
line), consisting of the stellar atmosphere having passed through 
the absorbing wind (dotted) and the thermal wind emission 
(dashed). Also included is the emission from a wind with pure 
free-free processes (long-dashed line) and the emission of the 
star without wind (dashed-dotted). 



same electron temperatures as the free-free Gaunt factors. The 
results are shown in the middle panel of Fig.|2] 

For A ^ 1 fim, the values of f(v, T) are «; 1 . The hydrogen 
free-bound processes, therefore, play no role in the IR and ra- 
dio regimes. At optical wavelengths however, i.e., for A < I fim, 
f(v, T) starts to grow steeply, with /(v, T) » 1 especially for 
electron temperatures below ~ 30 000 K, i.e., for values typically 
found in the winds of OB supergiants. This effect has been men- 
tioned akeady by Brussaard & van de Hulst (11962! ) and Kraus 
(12000) who pointed to the possible importance of the contribu- 
tion of free-bound processes to the total continuum emission. 

The resulting absorption coefficient, Ky, is shown in the bot- 
tom panel of Fig.|2]where we plot Ky/nl calculated with Eq. (fTST i 
for the different electron temperatures. Besides the well-known 
increase with increasing wavelength, the absorption coefficient 
also peaks at short wavelengths due to the growing influence of 
the free-bound processes. 

3. Results 

We fix the terminal wind velocity at Uoo = 1400 km s ' and the 
mass loss rate at M = 5 x 10"'' Moyr"'. In addition, we place the 
object to an arbitrary distance of 1 kpc. The resulting continuum 
emission is first calculated for the case of jS = 1.0 (Sect. 13. il l 
before we turn to the more interesting case of higher beta values 
in Sect.im 
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Fig. 4. Increase of the fluxes generated within wind zones with 
outer edge Ro^t, with respect to an infinitely large wind zone. The 
thermal wind emission in the optical and near-IR is generated 
completely in the vicinity of the stellar surface, i.e., within 5 
only. 

Table 1. Literature /3 values for OB supergiants. References are: 
M04 = Markova et al. (l2004] i: C06 = Crowther et al. ( [2006l l: 
K99 = Kudritzki et al. (1999"); E04 = Evans et al . (l2004] l: T05 = 
Trundle & Lennon (2005); T04 = Trundle et al. (l2004b . 



Galaxy 


Sp.Types 


/3 




Ref. 


Milky Way 


04- 


09.7 


0.7- 


1.25 


M04 


Milky Way 


09.5 


~B3 


1.2- 


3.0 


C06 


Milky Way 


BO- 


-B3 


1.0- 


3.0 


K99 


Magellanic Clouds 


08.5- 


-B0.5 


1.0- 


3.5 


E04 


Magellanic Clouds 


B0.5- 


-B2.5 


1.0- 


3.0 


T05 


Magellanic Clouds 


B0.5 


-B5 


1.0- 


3.0 


T04 



the wind. In Fig.|3] we also included the results for a wind with 
pure free-free emission. From a comparison of the total wind 
emission with the pure free-free wind emission, it is obvious that 
free-bound processes dominate the wind emission in the optical 
and near-IR wavelengths, i.e., for A ^ 2 fim. 

To understand where the optical wind continuum is gener- 
ated, we re-calculate the wind emission for different values of 
the outer edge of the wind zone. We then compared the resulting 
emission to the emission from an infinitely large wind, and we 
plotted the ratios versus frequency in Fig.|4] From this plot, it is 
evident that the wind emission in the optical and near-IR is gen- 
erated in the vicinity of the stellar surface, i.e., within 5 Most 
of the emission for A < lOjum (about 95 %) is even generated 
within the innermost 1.5 



3.1. A wind with 1.0 

For many hot star winds, /? lies typically in the range 0.8 ... 1 .0 
(e.g. Markova et al. 120041 Repolust et al. .2004) , and we use 
y6 = 1.0 for our first test calculation, which serves as a refer- 
ence model. The resulting spectrum is shown in Fig. [3] where 
we plot the total continuum emission of our test OB supergiant, 
consisting of the stellar atmosphere having passed through the 
absorbing wind and the thermal wind emission. For compari- 
son, we also included the emission of the star without wind. It 
is obvious, that the wind influences the total spectrum only for 
/I ^ 1 im\, while the optical spectrum remains uninfluenced. The 
stellar continuum only suffers in the IR and radio range, where 
for A 'Z lOjum the stellar emission is completely absorbed by 



3.2. Winds witli/3> 1.0 

In recent years, detailed investigations of the wind parameters 
of OB supergiants have revealed that for many objects the /3 pa- 
rameter in the velocity law given by Eq. (|2]) varies over a much 
larger range. Even values as high as 3.5 are reported. A list of 
/3 values of OB supergiants found in the literature is provided in 
Table[T] 

The main effect of a higher /3 value is the less steep increase 
in the wind velocity with distance from the stellar surface. This 
is shown in Fig.|5]where we plotted the velocity increase within 
5 Rt from the stellar surface in terms of the terminal velocity for 
different values of /3. According to the equation of mass conti- 
nuity (Eq. ([U), the wind density is proportional to v(ry \ Winds 
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0.5 




Fig. 5. Velocity increase in winds with different yS values. The 
higher the /?, the more slowly the wind is accelerated and the 
farther away from the star it reaches its terminal value. 
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Fig. 6. Density ratios close to the stellar surface for winds with 
y6 > 1 .0 with respect to the wind with /3 - 1.0. Due to the slower 
wind acceleration for higher the wind density remains much 
higher over a larger wind volume, resulting in a relative density 
enhancement close to the stellar surface, i.e., within 2 . . . 3 



with higher/? values consequently have a (much) higher density 
in the accelerating wind regions. 

We calculated the density distributions in winds with in- 
creasing /? and plotted the densities in terms of the density distri- 
bution for the wind with /3 = 1 .0 in Fig.|6] Close to the stellar sur- 
face the densities in the winds with yS > 1 are (much) higher than 
the reference value provided by the wind with yS = 1 . The higher 
the p, the stronger become these relative density enhancements 
due to (much) slower wind acceleration (see Fig.|5]l. These rela- 
tive density enhancements extend over the innermost ~ 2 . . . 3 /?* 
of the wind, while for larger distances, where the velocity distri- 
butions for the higher /3 values also start to approach the terminal 
velocity, the densities converge. 

In Sect. l3.1l we showed that for a wind with /? = 1.0 the op- 
tical continuum is generated within < 2Rt. This is exactly the 
range where /3 has its strongest influence on the wind density. If 
the optical continuum for winds with p > 1 .0 is generated within 
approximately the same region, we can expect that the resulting 
free-bound emission will be strongly influenced by the value of 

We determined the wind volume in which the optical flux is 
generated by comparing the flux generated within an infinitely 
large wind zone with the flux generated within winds with dif- 
ferent outer edges, /?out- The results for winds with different val- 
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Fig. 7. Flux ratio f'ff-fb(^co)//^if-fb(^out) as a function of outer 
wind edge R^ut for different values of The calculations shown 
are for A = 1 jjm and are also valid for all smaller wavelengths. 
As reference, we include the ratio for /?out = ^oo (solid line). For 
all values of the optical wind continuum emission is generated 
well within the plotted wind size of 5 R*. 
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Fig. 8. Fraction of stellar radiation escaping from the wind zone 
for different values of /3. While for yS = 1 .0 all the wind contin- 
uum emission for A < 1 yum passes unabsorbed through the wind, 
the situation is different for higher /3 values for which even the 
optical wind continuum suffers from wind absorption along its 
way through the wind. The higher/? the more stellar emission is 
absorbed also from the near- and mid-IR stellar spectrum. 



ues of /3 are shown as a function of the outer edge, /?out, in Fig.|2] 
These results are evaluated for A = 1 fim and are also valid for 
all smaller wavelengths (see Fig.|4|i. Obviously, for winds with 
/3 > 1.0, more than 90% of the optical continuum is generated 
within 1 .5 Rt, and more than 98% are produced within 2 R,. 

Before we turn to the calculation of the total continuum 
emission from winds with p > 1.0, we first investigate the in- 
fluence of yS on the stellar emission. Since the wind absorption 
coeflicient, /Cff_fb, is proportional to nl (see Eq. (fT2li). an increas- 
ing density leads to a (strongly) increasing optical depth. The 
stellar flux will therefore be more strongly absorbed in the wind 
with a high yS compared to a wind with /3 - I - This is illustrated 
in Fig. [8] Of course, the strongest depression of the stellar emis- 
sion occurs at IR wavelengths (A > I fim); with increasing/? also 
some part of the stellar flux at optical wavelengths, especially in 
the red part of the spectrum, is also absorbed. 

But while the increasing absorption coefficient results in a 
decrease of the stellar flux, it leads at the same time to a strong 
increase of the free-free and free-bound emission in the optical 
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Wavelength A [//m] 

Fig. 9. Comparison of total continuum spectra of the star plus 
wind system (solid lines) for different values of p. As in Fig.|3] 
a wind with pure free-free emission is included in each panel 
(long-dashed line) to emphasize the growing importance of the 
free-bound contributions for A < 2 /im. With increasing /3 the 
free-free and free-bound contributions (dashed) increase espe- 
cially in the near-IR and optical spectrum. At the same time, the 
stellar continuum (dotted) suffers from the increasing absorptiv- 
ity of the wind zone. For comparison, the pure stellar continuum 
is included (dashed-dotted). 



and IR part of the spectrum. This can be seen upon inspection of 
Fig. |9] where we plotted the total continuum emission for winds 
with different /3 values. The optical wind continuum increases 
with j6 leading to an enhanced total continuum in the IR and 
even at optical wavelengths. This increase, especially at optical 
wavelengths, is more evident in Fig.[TOl There we plot the flux 
ratio between the total continuum emission of the star plus wind 
system with /3 > 1 and the total continuum emission of the sys- 
tem with y6 = 1 as a function of wavelength. For y6 ^ 3.0, the total 
continuum exceeds the one fory6 = 1 in the optical by about 15% 
in the blue spectral range and by about 25-30% in the red spec- 
tral range, while in the near- and mid-IR it reaches even more 
than 200%. 

The contribution of the pure wind emission to the total con- 
tinuum production is finally shown in Fig.[TT] Its importance 
grows and shifts to lower wavelengths with increasing /3. While 
for /3 - 1.0 the contribution to the optical is less than 1%, it 




1.5 - 



Wavelength A [/^m] 



Fig. 10. Continuum flux ratio of the star plus wind system for 
different values of /3 with respect to yS = 1 .0. With increasing /3 
the increase in optical flux due to the increasing importance of 
the free-bound emission from the wind becomes visible. 




0.5 - 



Wavelength A [/inn] 

Fig. 11. Ratio of the wind emission with respect to the total con- 
tinuum emission for different values of /3. While for jS = 1 .0 the 
wind contribution to the total optical continuum is less than 1 %, 
its influence grows quickly with increasing /3, reaching values up 
to 30% in the red part of the optical spectrum for j0 > 3.0. 



grows to > 10 % for (i - 2.0 and to values exceeding even 30% 
foryS> 3.0. 

In summary, we can therefore state that even at optical wave- 
lengths the wind plays a non-negligible role for the continuum 
emission of OB supergiants having wind velocity distributions 
with/3 > 1.0. 

4. Discussion 

For our calculations we made use of several assumptions and 
simplifications like a constant wind temperature, the core-halo 
approximation, LTE even for the bound-free processes, and the 
neglect of electron scattering. These are severe restrictions to 
keep the model as simple as possible. Here, we want to discuss 
the influence and possible consequences of these assumptions 
and simplifications on the model results. In addition, we want 
to address shortly the topic of wind clumping and its expected 
influence on our results. 

4.1. The wind temperature 

For our model calculations of stellar wind emission (see 
Sect. l2.2l i we made two severe assumptions concerning the tem- 
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perature: (i) we assumed the wind to be isothermal, and (ii) we 
used a rather high temperature, i.e. the temperature at r - 
as the global electron temperature in the wind. Here, we want to 
discuss why these two assumptions are reasonable. 

Our choice of the electron temperature at r - means that 
we are calculating a lower limit of the wind absorption coeffi- 
cient, Ky, because it is proportional to r^'^^ and the Gaunt factors 
increase with decreasing temperature (see Fig.|2]i. Consequently, 
according to the Eqs. (|4]i and (|9]), the emission caused by free- 
free and free-bound processes, equally increases with decreasing 
electron temperature. Our results calculated for the high electron 
temperature thus underestimate the real wind continuum emis- 
sion. The error in wind emission can be estimated by looking 
at the real expected temperature distribution within the region 
where the optical wind continuum is generated. 

During our investigations we found that for each value of /3, 
95 % of the wind continuum at optical wavelengths is generated 
in the vicinity of the stellar surface, i.e. within 1.40 ... 1.65 Rf 
This is shown in the top panel of Fig. [121 The solid horizontal 
line in this plot indicates the 95 % level of the emission, and the 
vertical lines indicate, for each the distance at which this level 
is reached. 

To see how,in a more realistic wind model, the temperature 
would have changed over the wind region in which the optical 
continuum emission is generated, we finally calculate the tem- 
perature distribution according to Eq. (O, with the velocity dis- 
tribution as defined by Eq. (|2|. The resulting temperature distri- 
butions for different values of /3 are shown in the bottom panel of 
Fig-El An increase in /3 results in a lower wind velocity at the 
same distance from the star and, after Eq. (O, in a higher wind 
temperature. Therefore, the higher the /3, the hotter the wind re- 
mains at the same location. 

For winds with /3 > 2.0, the drop in wind temperature is 
found to be less than 5 % within the free-free and free-bound 
emission generation zone. Such a small decrease in temper- 
ature results in only a tiny and therefore negligible increase 
in wind emission. The assumption of an isothermal wind with 
Te - Ts(Rt) is therefore well justified. 

The situation is different for winds with /3 - 1.0. Here, 
the drop in wind temperature is found to be on the order of 
20 %, which corresponds to a decrease in electron temperature 
by about 5 000 K for our chosen model star. In such a case, the 
temperature distribution in the wind is not negligible. The de- 
crease in electron temperature results in an increase of the wind 
absorption coefficient of about a factor of two (see Fig.|2l), and 
consequently to a noticeably enhanced wind emission. However, 
even with an enhanced wind emission, the stellar spectrum still 
clearly dominates the total optical continuum. For OB super- 
giants with /3 - 1 .0, the assumption of an isothermal wind with 
Te = Te{Rt) is therefore still an acceptable approximation, as 
long as the mass loss rate of the star is not extremely high, as in 
the case of Luminous Blue Variables or B[e] supergiants. 

4.2. Core-halo approximation 

The construction of our star plus wind model is often referred 
to as the core-halo approximation, i.e., it is assumed that the 
stellar atmosphere is in hydrostatic equilibrium, and the wind is 
treated seperately with a density distribution following from the 
mass continuity equation. With such a treatment the transition 
between the atmosphere and the wind is not properly accounted 
for. 

However, the intention of our research was to investigate the 
influence of one single parameter, namely the steepness of the 
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Fig. 12. Top: As Fig.|2lbut for the innermost 2 only. The solid 
horizontal line indicates the 95 % level, and the vertical lines 
the distances where, for winds of different /3 values, this level is 
reached. Bottom: Electron temperature distribution within 2 
for different values of /3, normalized to the maximum electron 
temperature. For /3 > 1.0, the inner parts of the winds remain 
much hotter and almost at a constant value compared to the wind 
with/3 = 1.0. 



velocity increase, yS, keeping the stellar and the remaining wind 
parameters fixed. Deviations in density distribution of the star 
plus wind system introduced by the use of the core-halo approx- 
imation will consequently appear in all our discussed models 
and will influence all our results in the same way. But they will 
not significantly alter our conclusions, which are mainly based 
on the comparison of the spectra (i.e. the ratios) for stars with 
different /3 values. 

4.3. Non-LTE effects and electron scattering 

In our calculations, we assumed the wind to be in LTE. This 
is a justified assumption for the free-free processes since they 
are collisional. For the bound-free processes the situation is 
less clear Here, non-LTE effects might play an important role. 
According to our non-LTE model atmosphere calculations (e.g. 
Kubat l2003l l. departure coefficients can be (much) larger than 
unity in the outer parts of the stellar atmosphere (especially for 
the ground level) and, consequently, they lead to a higher bound- 
free opacity. 

We tested the influence of non-LTE effects on our results in 
a qualitative way by increasing artificially the bound-free opac- 
ity. This results in an increase in wind optical depth and, ac- 
cording to Eq. ( [Tol l, to an increased attenuation of the stellar 
emission, while the wind emission remains largely unaffected. 
Consequently, the flux of the total continuum emission decreases 
with increasing bound-free opacity, i.e., increasing influence of 
non-LTE effects. At the same time, the importance of the wind 
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emission with respect to the total continuum increases. The ne- 
glect of non-LTE effects thus results in a lower limit of the wind 
contribution. 

A similar conclusion can be drawn from the neglect of elec- 
tron scattering. Since electron scattering attenuates the stellar 
light passing through the wind zone, it equally acts as an addi- 
tional opacity source. The consequences for the total continuum 
emission are therefore the same: a decrease in total emission 
(over the wavelength ranges that are dominated by the stellar 
emission), and simultaneously an increase of the wind emission 
with respect to the total continuum. 

For our model star plus wind system, we checked the elec- 
tron scattering optical depth in radial direction. We find that for 
the chosen stellar and wind parameters and the different fi values, 
the electron scattering optical depth is always smaller than unity. 
However, compared to the free-bound opacity, it is not negligi- 
ble. This means, that we calculated indeed upper limits for the 
stellar emission resulting in lower limits for the wind emission 
with respect to the total continuum. A proper treatment of non- 
LTE effects and electron scattering can therefore only confirm 
our results: the importance of the wind contribution to the total 
optical continuum emission. 

4.4. Wind dumping 

In recent years, observational and theoretical investigations pro- 
vided evidence of wind clumping (see e.g. Hillier 120051) . One of 
the most striking results found from detailed spectroscopic anal- 
yses (e.g. Cr owther et al. [20021 Hillier et al. |2003l Bouret et al. 
120031 12005I I was that if a wind is clumped, the mass loss rates 
inferred from spectroscopy might be lower on average by a fac- 
tor of 3. In addition, Hillier et al. (I2003I I found that clumping 
starts close to the photosphere, i.e., within the region in which 
the influence of yS is strongest. 

The wind emission calculated in our study is not only deter- 
mined by the fi parameter of the velocity distribution, but also by 
the mass loss rate of the star that was kept constant during our 
analysis. Many of the authors who derived the high p values (see 
Table[T]i also claim that their mass loss rates are derived under 
the assumption of unclumped winds and might well be a factor 
of 3 lower Therefore, the question arises how our results might 
change if we account for wind clumping. 

This paper is not aimed to study wind clumping in detail, 
instead we refer to the recent paper by Kraus et al. (2008 ) who 
investigated in more detail the influence of wind clumping on 
the optical continuum emission of OB supergiants. Their results 
can be summarized as follows: wind clumping, introduced into 
the calculations, e.g., by the filling factor approach of Hiller et 
al. (120031 ). results in a slight decrease in wind emission at optical 
wavelengths. This decrease, however, was found to be less than 
the increase of the wind emission due to a high (3 value com- 
pared to a wind with p - 1.0. We can, therefore, conclude that 
in clumped winds with high p values the effects discussed in this 
paper are still present, but probably slightly weakened. 

5. Conclusions 

We investigated the influence of the thermal wind emission pro- 
duced by free-free and free-bound processes to the total contin- 
uum of normal OB supergiants. While for winds with a veloc- 
ity distribution following a y6-law with (3 typically in the range 
of 0.8 to 1.0, no influence of the wind at optical wavelengths 
is expected, the situation can be different when p exceeds 1.0. 



High p values reaching even 3.5 have recently been found for 
many OB supergiants. Our investigations therefore concentrated 
on such high beta values and their influence on the wind contin- 
uum emission especially at optical wavelengths. 

We found that the wind emission in the optical is generated 
within 2 iJ„ only. This region is exactly the region where p has 
its highest influence on the wind density structure. Since with 
increasing p the wind is accelerated much more slowly than for 
P — 1.0, the density close to the stellar surface is strongly en- 
hanced, leading to an enhanced production of especially free- 
bound emission at optical wavelengths. At the same time the 
stellar emission, which passes through the wind on its way to 
the observer, is absorbed. These effects increase with increasing 
p. The total continuum of OB supergiants, for which p values 
higher than 1.0 are found, can thus contain non-negligible con- 
tributions of wind emission even at optical wavelengths. 
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